Phase coherence in one dimensional superconductivity by power-law hopping 
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A central result in one-dimensional (ID) superconductivity is that even at zero temperature quan- 
tum fluctuations destroy phase coherence. Here we put forward a mechanism which can restore phase 
coherence: power-law hopping. We study a ID attractive-?/ Hubbard model with power-law hop- 
ping by Abelian bosonization and density-matrix renormalization group (DMRG) techniques. The 
parameter that controls the hopping decay acts as the effective, non-integer spatial dimensionality 
doff. We show analytically that for any dcd > 1 at zero temperature, power-law hopping suppresses 
fluctuations and induces phase coherence, namely, long-range superconducting order. A detailed 
DMRG analysis fully supports these flndings. These results are also of direct relevance to quantum 
magnetism as our model can be mapped onto a spin-chain with power-law decaying couplings, which 
can be studied experimentally by cold ion-trap techniques. 



PACS numbers: 74.78.Na, 74.40.-n, 75.10.Pq 

According to the Mermin-Wagner-Hohenberg theorem 
[TJ [5] quantum and thermal fluctuations in low dimen- 
sions prevent the spontaneous breaking of a continuous 
symmetry. A paradigmatic example is a one-dimensional 
(ID) superconductor (SC), where fluctuations of the SC 
order parameter result in quasi long-range order at zero 
temperature, i.e., the algebraic decay of the order pa- 
rameter correlation function [3]. By contrast supercon- 
ducting long-range order (LRO), equivalent to phase co- 
herence in this context, occurs if the correlation function 
does not decay even for arbitrarily large distances. 

Quantum phase-slip (QPS) fluctuations (i.e., processes 
where the amplitude of the order parameter temporarily 
vanish in some region along the wire, allowing its phase 
to slip by 2tt [3HS1) play an especially important role, 
as they are believed to destroy quasi-LRO and induce 
a finite resistivity at low temperatures observed in SC 
wires whose diameter is much smaller that the SC co- 
herence length [cf. Refs. [7| [H] and references therein]. 
Interestingly, recent theoretical works have shown the 
possibility to stabilize a ID SC through the suppression 
of QPS fluctuations. This is achieved by a weak cou- 
pling of the wire to a dissipative environment (such as 
a weakly disordered metallic thin fllm [5J (TU] , a metallic 
dirty lead [TT], a graphene sheet [H], or two dissipative 
electrodes attached to the ends of the SC wire [131 E]) 
that suppresses fluctuations, including QPS, and induce 
LRO [IS] . Experimentally, restoration of phase coherence 
has been recently observed in thin Zn [THl [T7] and Al [15] 
nanowires by increasing the coupling of the wire to dissi- 
pative electrodes. Therefore one of the main theoretical 
challenges in the fleld is to identify mechanisms that are 
capable to inhibit the proliferation of QPS and restore 
phase coherence in ID. The increase of the effective spa- 



tial dimensionality is an appealing choice. Interestingly, 
in the context of non-interacting ID weakly disordered 
systems [1^1 [20] , it is well known that power-law hopping 
cx 1/1* — j|" (with a > 1/2) effectively mimics the prop- 
erties of a system in dcs = 2/ (2a — 1) spatial dimensions 
with short-range hopping. This effect seems to be robust 
to the presence of interactions [H] . 

The achievement of phase coherence in low dimensional 
SC is also relevant for applications: from the miniaturiza- 
tion of the SC circuits to the enhancement of the critical 
temperature in SC nanostructures and thin films |22H25j . 

In this Letter we study the stabilization of a ID SC by 
power-law single-particle hopping. We focus our study 
on the ID attractive-L/ Hubbard model with algebraically 
decaying hopping tim oc t/\l~m\", where a is the param- 
eter controlling the decay. We study the quantum phases 
of the system at zero temperature by analytical (Abelian 
bosonization and a variational approach) and numeri- 
cal density-matrix renormalization group (DMRG) tech- 
niques. Our main result is the restoration of LRO at 
zero temperature for a < 3/2, corresponding to doff > 1- 
Algebraic coupling occurs in a variety of physical sys- 
tems, such as Josephson junction arrays [55], materials 
with strong dipolar interactions |27j . and atoms in cavi- 
ties realizing effectively quantum spin chains with long- 
range exchange interactions 28, 29 . In the latter, a spin- 
dependent optical dipole force applied to a cold atom gas 
makes possible to engineer power-law anti-ferromagnetic 
interactions with < a < 3 [55] [22] • As we show below, 
spin chains with power-law exchange can be mapped onto 
a ID SC with power-law hopping, so the ideas and tech- 
niques we introduce here are of direct relevance for these 
problems as well. 

Model- We study the L-site spin-1/2 ID Hubbard 
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model with attractive interaction U and power-law hop- 
ping, 
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where the fermionic annihilation operator ci^a destroys an 
electron at site / in spin state (T(=t,i) and hi^a = c\ ^ci^cr 
is the fermionic number operator. The (real) long-range 
hopping amplitude ti„i = t/\l — m|" connects sites I and 
m, /i is a uniform chemical potential corresponding to N 
particle per spin, and U controls the attractive interac- 
tion strength. For hopping restricted to nearest neigh- 
bors, solved exactly in [30,, SC quasi-LRO dominates 
over the competing charge- density wave correlations, ex- 
cept at half-filling where both correlations are compa- 
rable. On the other hand, purely imaginary power-law 
hopping in a ID repulsive-U Hubbard model at half fill- 
ing, investigated in [3T] for a = 1 (des = 2), induces a 
Mott metal-insulator transition at a finite value of the 
interaction strength. 

From now on we focus on the region \U\ ^ t 
where the local attractive interaction in Eq. ([!]) dom- 
inates. In this regime, unpaired electrons are ef- 
fectively forbidden at sufficiently low energies, and 
only Cooper pairs c|.^c| ^|0) are stable configura- 
tions. We therefore project out the singly-occupied 
sites at order t/|C/| with the unitary transformation 
n^g = e^He-'^, with S = -i {H+ - H^-) / \U\ and 

H^^ = -J2l^m,atlmi^-f^ls)cl^Cjnanma and H^T = 

-J2i^„i,a^i'mni^acl^c^„ (1 - h„ia)- The procedure is 
similar to the usual one employed to obtain the t-J model 
[32] . Here we mention the final result, 
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which effectively is a (long-range) variant of the well- 
known short-range Bose-Hubbard model with hard-core 
bosons [3]. Here hi = hi^^ + hi^^ is the total number op- 



erator at site I and A 



is the creation operator 



for a Cooper pair at site The term in square brackets 
in the above equation arises from second-order virtual 
processes in the hopping tim and contains the basic in- 
gredients leading to stabilization of the SC ground state 
driven by power-law hopping. Note that the coupling 
A|A„i minimizes the energy of the system by delocaliz- 
ing the Cooper pairs (thus favoring a more robust SC). 
By contrast the density-density interaction (hi — 1) h^ 
is strongly frustrated by power-law hopping. Therefore 



the competing charge density wave phase cannot be sta- 
bilized. The crucial sign difference between these two 
contributions is related to the assumption of purely real 
hoppings tijn- 

We now introduce the framework of the Abelian 
bosonization [5]. As a first step, we take the limit of van- 
ishing lattice parameter a — ?> in Eq. ^ and define the 
density hi /a — )■ p{x) and pair-creation Aj/a — )■ A^ (x) 
operators in the continuum. We next introduce the rep- 
resentation [3] 
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where (x) and (f) (x) are bosonic fields slowly varying on 
the scale of a. They satisfy the canonical commutation 
relations [\7(j) (x) , (y)] — inS [x ~ y). The field 9 [x) is 
physically related to the phase of the SC order parameter 



in the original system via (A(x)) — {c\ ^c^ 



-i9{x)\ 



while the field (j) (x) is related to slow Cooper pair density 
fluctuations Sp{x) ~ — V(/)(a;) /tt. 

The bosonic representation Q and Q allows to ex- 
press the low-energy part of Hamiltonian ([2| as 
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The first line of this equation is the Luttinger liquid 
model, where K is the dimensionless Luttinger param- 
eter controlling the asymptotic decay of the correlation 
function (e*''(^)e-*^('^')) - |x - x'\~^^^, and u is the ve- 
locity of the ID acoustic plasmons [5]. Physically, the 
product uK corresponds to the superfluid stiffness of the 
ID SC and K/u is the compressibility. The dimension- 
less coefficient g is a non-universal quantity measuring 
the strength of the power-law hopping term. The nu- 
merical values of K, u and g cannot be obtained from 
the bosonization procedure. Therefore we expect quali- 
tative rather than quantitative predictions. Finally, we 
note that in Eq. (l5| we have neglected higher harmon- 
ics Q^^p['t>('^)-'l'{'^^ arising from the non-local density- 
density interaction in Eq. ([2| , since the field (x) be- 
comes strongly fluctuating due to frustration. Its overall 
effect can be accounted by a renormalization of K. 

To make further progress, in what follows we employ 
the framework of the self-consistent harmonic approxi- 
mation (SCHA) [33 . This method consists in introduc- 
ing a Gaussian ansatz 5o = 2^ X]q5(7^ (^) ^q^q the 
Euclidean action of the system where q = (fc, — Wm) and 
LOra = inTm is the bosonic Matsubara frequencies at 
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temperature T [M]. The functions g^^ (q) are unknown 
variational parameters which must be chosen to minimize 
the variational free energy Fvar — Fq + T {S — Sq)q, with 
Fq the free energy associated to S'o, and S the action 
corresponding to Eq. ([s]). The notation (. . . )o stands for 
the average with respect to the trial action Sq. Minimiz- 
ing Fvar with respect to 50 {^) yields the self-consistent 
equation [3l[9l[TT] 
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where r = x ~ x' . In the regime 1/2 < a < 3/2, L ^■ 
00, T — >■ 0, an approximate solution of this self-consistent 
equation, asymptotically correct in the limit k 0, is 
given by the expression 

"^k^+vikr-'. (7) 
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The last term ~ r]\k\^°'~^ encodes the effect of power- 
law hopping in the system and is crucial in the rest of 
the analysis. Replacing ([t]) into ^ yields the equation 
for r? HH [35] 
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where G (?) = ^ 

k = ka, f = r/a, fj = 
Euler Gamma function 



dk 



■ cos I 



TT-qa^'"^" / uK and T (z) is the 
Jlj. Note that for fc ^- 0, the 
term ~ T^jfej^" ^ dominates over the term in Eq. 
for 1/2 < a < 3/2. This is the key ingredient for 
the restoration of phase coherence in this region. For 
a < 1/2, the integral in Eq. (|6| becomes divergent due 
to the slowly decaying hopping, and consequently expres- 
sion ([7| is no longer a valid solution. In this regime, 
fluctuations around the ordered groundstate with fc = 
are fully suppressed and mean-field theory becomes ex- 
act. For a > 3/2 the term k^ dominates for fc -> 0, 
and the model can be effectively mapped onto a ID SC 
with renormalized short-range couplings [TJ 1301 • There- 
fore, ttc = 3/2 is the critical value that separates the 
regimes of quasi-LRO from robust superconducting LRO. 
In order to further characterize the role of power-law hop- 
ping in ID SC, we now compute the average of the order 
parameter (e*''(^))o = e-<'''(°)>''/2 = exp [-G (0) /2]. In 
stark contrast with the short range case |3J |3D] the aver- 
age SC order parameter is finite for any 1/2 < a < 3/2. 
Moreover, the equal-time pair correlation function be- 
comes C(r) = (^e'^{r)ed,it~^oio)^^ ^ exp [G (r) - G (0)], 

which in the limit r — > cx) tends to a constant C (r) « 
g-G(o) [1 _^ _^ Q ^^/^3-2a)j ^ ^ji-i^ ^ > 0. Di- 

rect comparison of (6'^ (0))o in our case and in a short- 
ranged doff-dimensional results in an explicit expression 
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FIG. 1: (Color online) Analytic results for C (r/a) as obtained 
with the SCHA for different values of a. In order to compare 
with DMRG results, K = 1.8, poa = 2N/L = 0.1803. For a = 
0.7 (continuous black line), the onset of SC LRO is evident 
in the emergence of a plateau from r/a « 10 (^ = 1.5), while 
for Q = 1.3 (inset) the characteristic length of the plateau is 
^ — 1.6 X 10^. Instead, quasi-LRO is observed for r < ^, with 
effective Luttinger parameter Kaft ~ 3.7. 



of doff as a function of a, dcs — 2/ {2a — 1). Therefore, 
for doff > 1, corresponding to 1/2 < a < 3/2, LRO and 
phase coherence are restored (see Fig. [T]). This is the 
main result of this Letter. Finally we estimate the min- 
imum length scale ^ at zero temperature necessary to 
observe LRO by equating the contributions |fc|^" ^ and 
P in Eq. (l7|. 
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where it is assumed that {L,r} ^ ^ [cf. Figjl]]. 
Numerical results.- We now test the results from the 
bosonization approach above by DMRG j571 - [5^ . Power- 
law hopping is a challenge for many-body numerical sim- 
ulations as finite size effects become much more impor- 
tant. In the context of DMRG also the number of basis 
states that must be kept increases dramatically with re- 
spect to short-range models. Moreover as the critical 
value ac = 3/2 is approached, the crossover length scale 
^ diverges [cf. Eq.(|9|], and DMRG is unable to reach the 
LRO region. With these limitations in mind, we compute 
the spatial average of the pair correlation function 



L-lo-r 
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using the DMRG, where (...) stands for the average in 
the ground state of the microscopic model Eq. ([!]), and 
Iq is the number of sites at the end of the chain which 
are eliminated in order to minimize finite-size effects. In 
the limit \U\ /t ^ 1, where amplitude fluctuations of the 
SC order parameter are negligible, the correlation C{r) 
[normalized by G(0)] should compare to the analytical 
results of Fig. [l] In Fig. [2] we plot C (r) as a function 
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FIG. 2: (Color online) C(r) Eq. ([To|, computed by DMRG, 
for L = 2S3,N = 21, U/t = -20, = L/4: and different a's. 
In agreement with the bosonization results, LRO is clearly 
observed up to a ~ 1. For larger a, the crossover length 
[cf. Eq. [9] to observe LRO is larger than the maximum 
size accessible by DMRG techniques. The standard deviation 
(obtained by taking the spatial average) is shown only for 
Of = 0.5. A finite size scaling analysis with almost constant 
N/L of C(r) (inset) in the region a ~ 1 indicates that the 
decay with the distance becomes slower as we increase the 
system size. It is therefore plausible that C (r) will eventually 
flatten out as the L — >■ 00 is approached. 



of r for different values of a and \U\ /t = 20. The nu- 
merical results are consistent with the previous analytical 
calculation. LRO seems to occur for a<l.l though as 
a increases a very slow downward trend is observed. In 
order to further clarify whether the very slow decay ob- 
served for a ~ 1 disappears in the thermodynamic limit 
we carry out a finite size scaling analysis. The results, 
depicted in Fig. [2] (inset), show that the decay becomes 
slower as the system size increases. This is a strong in- 
dication that true LRO also occurs for 1 < a < 3/2. 
Qualitatively similar results are obtained for other filling 
factors and interaction strength provided that \U\/t ^ 1, 
in order to avoid that £, > L. For a > 1.1, also in agree- 
ment with the analytical calculation, only quasi-LRO is 
observed for the sizes accessible by DMRG techniques. 
Much larger sizes L are needed to observe LRO. In con- 
clusion, DMRG results indicate that for \U\/t 3> 1 and 
a < 1-1 LRO is restored by power-law hopping. 

Implications for transport- While a detailed study 
of the transport properties is beyond the scope of the 
present Letter, it is clear that the last term in Eq. ([5| 
will have strong effects on the QPS/anti-QPS interac- 
tion, which determines the resistivity of the ID super- 
fluid [3 HO]- In our case, we expect the attractive inter- 
action between a QPS located at = (xi,uTi), and an 
anti-QPS located at rj = {xj,UTj) to grow faster than 
In (|ri — r^l /a), the typical interaction in a short-ranged 
ID SG [3110]. Therefore, in our case the QPS/anti-QPS 
pair should be tightly bound and consequently we expect 
the resistivity of the system to approach to zero much 
faster than the short-ranged case as T — )■ and L 00. 



A more detailed study of the tranport properties of the 
wire is currently in progress '41]. 

Application to quantum spin chains.- These findings 
not only provide novel insights into the nature of the su- 
perconducting state in low spatial dimensions, but are 
also of interest to quantum magnetism. Using the An- 
derson pseudo-spin representation [12] ni — > S'f -I- 1/2, 
Aj — !• , Eq. ([2]) can be mapped onto a spin-1/2 XX Z 
chain with an effective Zeeman field along the z-axis, and 
long-range antiferromagnetic (ferromagnetic) couplings 
along the z-axis (xy-plane) . In this form, we can see 
immediately that the long-range nature of the exchange 
couplings will induce frustration along the z-axis, but 
will favor ferromagnetic LRO in the xy-plane. As men- 
tioned above, it recently became experimentally feasible 
to engineer a broad range < a < 3 of power-law decay- 
ing interactions in cold trapped-ion systems that mimic 
quantum spin chains [281 129j . Our results are therefore 
relevant in the description of the rich phase digram of 
these systems. 

In conclusion, we have investigated the ID attractive- 
U Hubbard model with an algebraic hopping by means 
of Abelian bosonization and DMRG techniques. Results 
from both approaches are consistent: at T = 0, true 
LRO is recovered for a < 3/2, corresponding to an effec- 
tive dimensionality dcs > 1. However, for a > 1.1 LRO 
cannot be observed numerically because it requires sys- 
tem sizes far larger than those accessible by DMRG. The 
robustness of superconductivity in des > 1 paves also the 
way to boost superconductivity by shell effects and 
other coherence effects important in low-dimensional and 
nanoscale SGs. Our results are of interest in other prob- 
lems beyond superconductivity, especially in ID quan- 
tum magnetism, where the resulting phase diagram can 
be investigated experimentally in trapped-ion systems. 
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